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Abstract 

We calculate the nonperturbative J/^ — N and \&' — N cross sections with the 
model of the stochastic vacuum which has been succesfully applied in many 
high energy reactions. We also give a quantitative discussion of pre-resonance 
formation and medium effects. 
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During the next year the first data on heavy ion collisions at high energy (y/s = 200 GeV 
per nucleon pair) will be available at RHIC. As it is well known, one of the main goals of 
this machine is to find and study the plasma of quarks and gluons (QGP) The search 

for this state of matter has started in the early eighties and since then has been subject of 
intense debate |],|2]]. Recent results on Pb-Pb collisions, taken at lower energies at CERN- 
SPS, have attracted a great attention and increased the hope that a new phase of nuclear 
matter is "just around the corner" . 

The signature of QGP formation has been and still remains a theoretical and experi- 
mental challenge. Indeed there is so far no "crucial test" able to disentangle the possible 
new phase from the dense hadronic background. Among the proposed signatures the most 
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interesting is the suppression of J/ty jj. This suppression was observed experimentally by 
the NA38 collaboration in 1987 in collisions with light ions and also and more dramatically 
by the NA50 collaboration in 1995-1996 in Pb-Pb collisions [g||. 

Whereas the old data could be reasonably well explained by a "conventional" approach 
the new Pb-Pb results of the NA50 collaboration created a big controversy |J . The entire 
set of the J/^> data from pA and AB collisions available before the advent of the Pb beam 
at CERN SPS has been found to be consistent with the nuclear absorption model. However, 
in the case of the new Pb-Pb data, the density of secondaries, which (together with primary 
nucleons flying around) are presumably responsible for the charmonium absorption, is so 
high that the hadronic system in question is hardly in a hadronic phase. Nevertheless, a 
conventional treatment of the problem is not yet discarded |7j] . 

Reliable values for the charmonium nucleon cross sections are of crucial importance 
in the present context. One needs to know the cross section aj/^ in order to predict a 
nuclear suppression of J/\I/ without assuming a so called "deconfining regime". Estimates 
using perturbative QCD give values which are too small to explain the observed absorption 
conventionally, but they are certainly not reliable for that genuine nonperturbative problem. 
A nonperturbative estimate may be tried by applying vector dominance to J/^ and ^' 
photoproduction. In this way a cross section of (Jji^ — 1.3 mb for ^/s ~ 10 GeV and 



oy/erjA/, — 0.8 has been obtained [§,0. In ref. |T0j strong arguments against the vector 
dominance with only few intermediate vector mesons were put forward even in the case of 
the production of light vector mesons and which apply a fortiori to the production of heavy 
vector mesons. The instability of the vector dominance model can be seen from a more 
refined multichannel analysis || where a value <Jj/^p — 3 — 4 mb has been obtained. Even 
this value is too small in order to explain the absorption in p-A collisions which is of the 



order aff s ~ 7.3 mb |]TT],|T^ 



These hadron-hadron cross sections involve nonperturbative aspects of QCD dynamics 
and therefore require a nonperturbative model to be calculated. In a recent letter |T3J the 
nonperturbative QCD contribution to the charmonium-nucleon cross section was evaluated 
by using an interpolation formula for the dependence of the cross section on the transverse 
size of a quark-gluon configuration. In this work we calculate the J/\l/ and - nucleon cross 
sections in a specific nonperturbative model of QCD: the model of the stochastic vacuum 
(MSV) fl4}-|T7|. It has been applied to a large number of hadronic and photoproduction 
processes (including photoproduction of J /^) with remarkably good success. Its application 
to J/\l/ and ^' nucleon scatttering is straightforward. We also investigate the influence of the 
nuclear matter and arrive at rather stringent limits for the cross sections in an environment 
different from the vacuum and where the properties of the medium are reflected in a shift 
of the J/\l> mass |TJ]. 

The basis of the MSV is the calculation of the scattering amplitude of two colourless 
dipoles ||T0| , |T7| based on a semiclassical treatment developped by Nachtmann For de- 
tails we refer to the literature and show here only some intermediate steps necessary for the 
understanding of the text. The dipole-dipole scattering amplitude is expressed as the ex- 
pectation value of two Wegner- Wilson loops with lightlike sides and transversal extensions 
fa and f t 2 respectively. This leads to a profile function J(b, fti, r^) from which hadron- 
hadron scattering amplitudes are obtained by integrating over different dipole sizes with the 
transversal densities of the hadrons as weight functions according to 
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tot 



d 2 bd 2 r tl d 2 r t2 pj/q,(f tl ) p N (r t2 ) J(b,r tl ,r t2 ) . (1) 



Here pj/y(r tl ) and pNift2) are the transverse densities of the J/^f and nucleon respectively. 

The basic ingredient of the model is the gauge invariant correlator of two gluon field 
strength tensors. The latter is characterized by two constants: the value at zero distance, 
the gluon condensate < g 2 FF >, and the correlation length a. We take these values from 
previous applications of the model ( and literature quoted there): 

< g 2 FF >= 2.49GeV 4 a = 0.346fm . (2) 

The wave functions of the proton have been determined from proton-proton and proton- 
antiproton scattering respectively. It turns out that the best description for the nucleon 
transverse density is given by that of a quark diquark system with transversal distance r t 
and density: 

PN (r t ) = \^ p (f t )\ 2 = ^e^ . (3) 

The value of the extension parameter, S p = 0.739 fm, obtained from proton-proton scatter- 
ing agrees very well with that obtained from the electromagnetic form factor in a similar 
treatment |p0[| . 

For the wave function of the J/\l/ we used two approaches: 
1) A numerical solution of the Schroedinger equation with the standard Cornell potential 



21 : 



V = -\^ + ar. (4) 
3 r 

2) A Gaussian wave function determined by the electromagnetic decay width of the J/ty 
which has been used in a previous investigation of J/\I/ photoproduction fT0|j . 

For the no analysis of photoproduction in the model has been made so we use only 
the solution of the Schroedinger equation. 

The linear potential can be calculated in the model of the stochastic vacuum which yields 
the string tension: 

a = — ^ < g 2 FF > a 2 = 0.179 GeV 2 , (5) 

817T 



where the parameter k has been detemined in lattice calculations to be k = 0.8 [f2~!2] 



The other parameters, the charmed (constituent) mass and the (frozen) strong coupling 
can be adjusted in order to give the correct J/\I/ and mass difference and the J/ip decay 
width 

m c = 1.7GeV a s = 0.39 . (6) 



We also use the standard Cornell model parameters ||21|| : 

a s = 0.39, a = 0.183 GeV 2 and m c = 1.84 GeV. (7) 
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From the numerical solution VKI^I) °f the Schroedinger equation the transversal density 
is projected: 



pj/*{r tj 



ip(yjr t 2 + r$) 



dr 



3 j 



(8) 



where r t is the J/\I/ transversal radius. 

Given the values of a s , a and m c we solve the non-relativistic Schroedinger equation 
numerically, obtain the wave function, compute the transverse wave function and plugg it 



into the MSV calculation fl7|l. The results are shown in Table I. In this table \[<r 2 > is 



the root of the mean square distance of quark and antiquark and y < rf > is the root of 
the mean square transversal distance of quark and antiquark. Wave function A) is the one 
obtained with the parameters given by Eqs. (||) and (||). Wave function B) corresponds to 
the standard Cornell model parameters [^T|] Eq. (|7]). 



Wave function 


V< r 2 > fm 


sj '< r\ > fm 


<?tot [mb] 


J/V(1S) 








A 


0.393 


0.321 


4.48 


B 


0.375 


0.306 


4.06 


C 






4.69 


#(2£) 








A: 


0.788 


0.640 


17.9 



TABLE I J/*S> — N and ty' — N cross section. A and B: numerical solution of the Schroedinger 

equation with parameters in Eqs. (^) and (^) respectively. C: Cross section obtained by the 
weighted average of the longitudinally and transversely polarized J/^> wave functions of ref. [10|. 



In ref. [ 10 1 a gaussian ansatz was made to construct vector meson wave functions that 



describe well the electromagnetic decay of the vector meson and photo and electroproduction 
cross sections. In Table I, wave function C) gives the result for the J — N cross section 
obtained with the weighted average of the longitudinally and transversely polarized J/ty 

0.327 fm and 0.466 fm. 



wave functions from [JTOj] with transversal sizes y < rf > 

Averaging over our results for different wave functions, our final result for the J/ty — N 
cross section is 



4.4±0.6mb 



(9) 



The error is an estimate of uncertainties coming from the wave function and the model. 
The only other nonperturbative calculation of the J /\& — N cross section that we are aware 
of was done in ref. |H| and the obtained cross section was aj/^ = 3.6 mb, in a fair agreement 



with our result and with recent analysis of J/\l/ photoproduction data M. For our cross 



section is also of the same order as the value obtained in |13|| : oy = 20.0 mb. 

Since one of the possible explanations of the observed J/\l/ suppression is based on 
the pre-resonance absorption model |Tl| we present numerical calculations of the nucleon - 
pre-resonant charmonium state cross section, a^. In the pre-resonance absorption model, 
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the pre-resonant charmonium state is either interpreted as a color-octet, (cc)g, and a 
gluon in the hybrid (cc)g — g state, or as a coherent J/\& — \P' mixture. We use a gaus- 
sian transverse wave function, as in Eq. ([!]), to represent a state with transversal radius 

\J < r 2 > ~ 0.82 V < > = V^SV (^0 * s ^ ne P r e-resonance extension parameter analogous 
to S p ). With the knowledge of the wave functions and transformation properties of the con- 
stituents we can compute the total cross section given by the MSV. The resulting nucleon - 
pre-resonant charmonium state cross section will be different if the pre-resonant charmonium 
state consists of entities in the adjoint representation (as (cc) 8 — g) or in the fundamental 

representation (as a J/$ - $' mixture), the relation being <7 a dj mt = N z f i ^fundamental ; with 
Nc = 3. In Table II we show the results for these two possibilities and different values of 
the transverse radius. 



(fm) 


(mb) 


a {ce)8-g 

(mb) 


0.20 


1.79 


4.02 


0.25 


2.76 


6.21 


0.30 


3.96 


8.91 


0.35 


5.30 


11.92 


0.40 


6.81 


15.32 


0.45 


8.50 


19.12 


0.50 


10.28 


23.13 



TABLE II The cross section charmonium-nucleon for gaussian wave-functions and different 
values of the transverse radius (J < r 2 >) of the cc in a singlet state (first row) or in a hybrid 

(cc)s — g state (second row). 



From our results we can see that a cross-section ~ 6 — 7 mb, needed to explain the 
J/\P and suppression in p-A collisions in the pre-resonance absorption model |TT|,|7|, is 
consistent with a pre-resonant charmonium state of size ~ 0.50 — 0.55 fm if it is a J/ty — ^' 
mixture or ~ 0.30 — 0.35 fm for a (cc) 8 — g state. 

So far the calculations were done with the vacuum values of the correlation length and 
gluon condensate generally used in the MSV [17] . However, since the interaction between the 
charmonium and the nucleon occurs in a hadronic medium, these values may change. Indeed, 
lattice calculations |2"2"| , |2T| show that both the correlation length and the gluon condensate 
tend to decrease in a dense (or hot) medium. The reduction of the string tension, a, leads to 
two competing effects, which can be quantitatively compared in the MSV. On one hand the 
cross section tends to decrease strongly when the gluon condensate or the correlation length 
decrease. On the other hand, when the string tension is reduced the c — c state becomes less 
confined and will have a larger radius, which, in turn, would lead to a larger cross section 
for interactions with the nucleons in the medium. It is of major interest to determine which 
of these effects is dominant. 

The dependence of the total cross section, Eq. ([!]), on the extension parameters S p and 
S^p is quite well parametrized [Q as: 
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a m oc< g FF > z a iU ( ^ pi 



Q \ 1-5 / C s 2 



(10) 



In the MSV the string tension, cr, is related to the gluon condensate and to the correlation 
length through Eq. (^j). Therefore, the dependence of the cross section on the string tension 
and correlation length is approximately given by: 



(f f (t ' 



(id 



In a rough approximation the hadron radii can be estimated, using the Ritz variational 
principle, to be: 



(12) 



and thus we finally obtain the following three possibilities to express the cross section as 
a function of the string tension, a, the correlation length, a, and the gluon condensate, 

< g 2 FF >: 

a^ N oc \ a 25 / 12 < g 2 FF >- 5 / 4 (13) 
{ < g 2 FF > 5//6 a 25//6 

From the equation above we see that the final effect of the medium is a reduction in the 
cross section. We can also see that a 10% variation in the parameters lead to large variations 
on the cross sections. Using the values of the correlation length and the gluon condensate 
reduced by 10%: a = 0.31 fm , (g 2 FF) = 2.25 GeV 4 , we obtain a 40% reduction in the 
cross sections. Taking this reduction into account the cross sections obtained in this work 
are smaller than the ones needed (both in Refs. |IT| and J/J) to explain experimental data. 
However, since in our model the nucleon - pre-resonant charmonium state cross section is 
much bigger than the J/^J — N cross section if the pre-resonance charmonium state is a 
hybrid (cc) 8 — g state, the reduction in the cross sections due to medium effects favors the 
pre-resonant model for the hadronic explanation of the observed J/^f suppression. 

In order to get more precise results we have varied only one the parameters a and 

< g 2 FF > and kept the other fixed. This was done in such a way as to decrease the string 
tension according to equation (5) to the values given in Table III, first row. The numerically 
evaluated values for the mass of the J/^f (serving as a physical measure of the change) and 
the cross sections are given in rows 2 to 4, using Gaussian wave functions determined by the 
variational principle. 



string tension [GeV 2 ] 


AE [MeV] 


<rtot [mb] 


crtot [mb] 


(GeV 2 ) 


MeV 


(mb) 


(mb) 






a const. 


< g 2 FF > const. 







4.48 


4.48 


0.9(7o 


-33 


3.86 


3.38 


0.75(7 


- 83 


3.27 


2.2 


0.50(7o 


-174 


2.0 


0.83 


0.25ao 


-275 


0.91 


0.13 
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TABLE III J/^-Nucleon total cross sections as function of the string tension with either the 
correlation length or the gluon condensate kept constant. AE is the mass decrease of the J/^f 

due to the change of string tension. 



To summarize, we calculated the nonperturbative J /\& — iV and ty' — N cross sections 
with the MSV. The basic ingredient of the model is the gauge invariant correlator of two 
gluon field strength tensors which is characterized by two constants: the gluon condensate 
and the gluon field correlation length. Using for these quantities values fixed in previous 
applications and using well accepted charmonium wave functions we obtain uj^ n ~ 4 mb 
and ay jv ~ 18 mb. An interesting prediction of the MSV is the strong depedence on the 
parameters of the QCD vacuum which will most likely lead to a drastic reduction of them 
at higher temperatures and perhaps also at higher densities. 
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